Studying Plant Specialized Metabolites Using Computational Metabolomics Strategies.

Plant specialized metabolites have diversified vastly over the course of plant evolution, and they are considered key players in complex interactions between plants and their environment. The chemical diversity of these metabolites has been widely explored and utilized in agriculture and crop enhancement, the food industry, and drug development, among other areas. However, the immensity of the plant metabolome can make its exploration challenging. Here we describe a protocol for exploring plant specialized metabolites that combines high-resolution mass spectrometry and computational metabolomics strategies, including molecular networking, identification of structural motifs, as well as prediction of chemical structures and metabolite classes.

POMBOX: A Fission Yeast Cloning Toolkit for Molecular and Synthetic Biology.

The fission yeast Schizosaccharomyces pombe is a popular model organism in molecular biology and cell physiology. With its ease of genetic manipulation and growth, supported by in-depth functional annotations in the PomBase database and genome-wide metabolic models,S. pombe is an attractive option for synthetic biology applications. However,S. pombe currently lacks modular tools for generating genetic circuits with more than 1 transcriptional unit. We developed a toolkit to address this gap. Adapted from the MoClo-YTK plasmid kit for Saccharomyces cerevisiae and using the same modular cloning grammar, our POMBOX toolkit is designed to facilitate fast, efficient, and modular construction of genetic circuits inS. pombe. It allows for interoperability when working with DNA sequences that are functional in bothS. cerevisiae and S. pombe (e.g., protein tags, antibiotic resistance cassettes, and coding sequences). Moreover, POMBOX enables the modular assembly of multigene pathways and increases the possible pathway length from 6 to 12 transcriptional units. We also adapted the stable integration vector homology arms to Golden Gate assembly and tested the genomic integration success rates depending on different sequence sizes, from 4 to 24 kb. We included 14 S. pombe promoters that we characterized using two fluorescent proteins, in both minimally defined (EMM2─Edinburgh minimal media) and complex (YES─yeast extract with supplements) media. Then, we examined the efficacy of 6 S. cerevisiae and 6 synthetic terminators in S. pombe. Finally, we used the POMBOX kit for a synthetic biology application in metabolic engineering and expressed plant enzymes in S. pombe to produce specialized metabolite precursors, namely, methylxanthine, amorpha-4,11-diene, and cinnamic acid from the purine, mevalonate, and aromatic amino acid pathways.

Machine Learning-Guided Protein Engineering.

Recent progress in engineering highly promising biocatalysts has increasingly involved machine learning methods. These methods leverage existing experimental and simulation data to aid in the discovery and annotation of promising enzymes, as well as in suggesting beneficial mutations for improving known targets. The field of machine learning for protein engineering is gathering steam, driven by recent success stories and notable progress in other areas. It already encompasses ambitious tasks such as understanding and predicting protein structure and function, catalytic efficiency, enantioselectivity, protein dynamics, stability, solubility, aggregation, and more. Nonetheless, the field is still evolving, with many challenges to overcome and questions to address. In this Perspective, we provide an overview of ongoing trends in this domain, highlight recent case studies, and examine the current limitations of machine learning-based methods. We emphasize the crucial importance of thorough experimental validation of emerging models before their use for rational protein design. We present our opinions on the fundamental problems and outline the potential directions for future research.

transXpress: a Snakemake pipeline for streamlined de novo transcriptome assembly and annotation.

BACKGROUND: RNA-seq followed by de novo transcriptome assembly has been a transformative technique in biological research of non-model organisms, but the computational processing of RNA-seq data entails many different software tools. The complexity of these de novo transcriptomics workflows therefore presents a major barrier for researchers to adopt best-practice methods and up-to-date versions of software. RESULTS: Here we present a streamlined and universal de novo transcriptome assembly and annotation pipeline, transXpress, implemented in Snakemake. transXpress supports two popular assembly programs, Trinity and rnaSPAdes, and allows parallel execution on heterogeneous cluster computing hardware. CONCLUSIONS: transXpress simplifies the use of best-practice methods and up-to-date software for de novo transcriptome assembly, and produces standardized output files that can be mined using SequenceServer to facilitate rapid discovery of new genes and proteins in non-model organisms.

Guiding the choice of informatics software and tools for lipidomics research applications.

Progress in mass spectrometry lipidomics has led to a rapid proliferation of studies across biology and biomedicine. These generate extremely large raw datasets requiring sophisticated solutions to support automated data processing. To address this, numerous software tools have been developed and tailored for specific tasks. However, for researchers, deciding which approach best suits their application relies on ad hoc testing, which is inefficient and time consuming. Here we first review the data processing pipeline, summarizing the scope of available tools. Next, to support researchers, LIPID MAPS provides an interactive online portal listing open-access tools with a graphical user interface. This guides users towards appropriate solutions within major areas in data processing, including (1) lipid-oriented databases, (2) mass spectrometry data repositories, (3) analysis of targeted lipidomics datasets, (4) lipid identification and (5) quantification from untargeted lipidomics datasets, (6) statistical analysis and visualization, and (7) data integration solutions. Detailed descriptions of functions and requirements are provided to guide customized data analysis workflows.

Ion identity molecular networking for mass spectrometry-based metabolomics in the GNPS environment.

Molecular networking connects mass spectra of molecules based on the similarity of their fragmentation patterns. However, during ionization, molecules commonly form multiple ion species with different fragmentation behavior. As a result, the fragmentation spectra of these ion species often remain unconnected in tandem mass spectrometry-based molecular networks, leading to redundant and disconnected sub-networks of the same compound classes. To overcome this bottleneck, we develop Ion Identity Molecular Networking (IIMN) that integrates chromatographic peak shape correlation analysis into molecular networks to connect and collapse different ion species of the same molecule. The new feature relationships improve network connectivity for structurally related molecules, can be used to reveal unknown ion-ligand complexes, enhance annotation within molecular networks, and facilitate the expansion of spectral reference libraries. IIMN is integrated into various open source feature finding tools and the GNPS environment. Moreover, IIMN-based spectral libraries with a broad coverage of ion species are publicly available.

Feature-based molecular networking in the GNPS analysis environment.

Molecular networking has become a key method to visualize and annotate the chemical space in non-targeted mass spectrometry data. We present feature-based molecular networking (FBMN) as an analysis method in the Global Natural Products Social Molecular Networking (GNPS) infrastructure that builds on chromatographic feature detection and alignment tools. FBMN enables quantitative analysis and resolution of isomers, including from ion mobility spectrometry.

The chloroalkaloid (-)-acutumine is biosynthesized via a Fe(II)- and 2-oxoglutarate-dependent halogenase in Menispermaceae plants.

Plant halogenated natural products are rare and harbor various interesting bioactivities, yet the biochemical basis for the involved halogenation chemistry is unknown. While a handful of Fe(II)- and 2-oxoglutarate-dependent halogenases (2ODHs) have been found to catalyze regioselective halogenation of unactivated C-H bonds in bacteria, they remain uncharacterized in the plant kingdom. Here, we report the discovery of dechloroacutumine halogenase (DAH) from Menispermaceae plants known to produce the tetracyclic chloroalkaloid (-)-acutumine. DAH is a 2ODH of plant origin and catalyzes the terminal chlorination step in the biosynthesis of (-)-acutumine. Phylogenetic analyses reveal that DAH evolved independently in Menispermaceae plants and in bacteria, illustrating an exemplary case of parallel evolution in specialized metabolism across domains of life. We show that at the presence of azide anion, DAH also exhibits promiscuous azidation activity against dechloroacutumine. This study opens avenues for expanding plant chemodiversity through halogenation and azidation biochemistry.

Metabolomics Data Preprocessing Using ADAP and MZmine 2.

The informatics pipeline for making sense of untargeted LC-MS or GC-MS data starts with preprocessing the raw data. Results from data preprocessing undergo statistical analysis and subsequently mapped to metabolic pathways for placing untargeted metabolomics data in the biological context. ADAP is a suite of computational algorithms that has been developed specifically for preprocessing LC-MS and GC-MS data. It consists of two separate computational workflows that extract compound-relevant information from raw LC-MS and GC-MS data, respectively. Computational steps include construction of extracted ion chromatograms, detection of chromatographic peaks, spectral deconvolution, and alignment. The two workflows have been incorporated into the cross-platform and graphical MZmine 2 framework and ADAP-specific graphical user interfaces have been developed for using ADAP with ease. This chapter summarizes the algorithmic principles underlying key steps in the two workflows and illustrates how to apply ADAP to preprocess LC-MS and GC-MS data.

PBS3 and EPS1 Complete Salicylic Acid Biosynthesis from Isochorismate in Arabidopsis.

Salicylic acid (SA) is an important phytohormone mediating both local and systemic defense responses in plants. Despite over half a century of research, how plants biosynthesize SA remains unresolved. In Arabidopsis, a major part of SA is derived from isochorismate, a key intermediate produced by the isochorismate synthase, which is reminiscent of SA biosynthesis in bacteria. Whereas bacteria employ an isochorismate pyruvate lyase (IPL) that catalyzes the turnover of isochorismate to pyruvate and SA, plants do not contain an IPL ortholog and generate SA from isochorismate through an unknown mechanism. Combining genetic and biochemical approaches, we delineated the SA biosynthetic pathway downstream of isochorismate in Arabidopsis. We found that PBS3, a GH3 acyl adenylase-family enzyme important for SA accumulation, catalyzes ATP- and Mg2+-dependent conjugation of L-glutamate primarily to the 8-carboxyl of isochorismate and yields the key SA biosynthetic intermediate, isochorismoyl-glutamate A. Moreover, we discovered that EPS1, a BAHD acyltransferase-family protein with a previously implicated role in SA accumulation upon pathogen attack, harbors a noncanonical active site and an unprecedented isochorismoyl-glutamate A pyruvoyl-glutamate lyase activity that produces SA from the isochorismoyl-glutamate A substrate. Together, PBS3 and EPS1 form a two-step metabolic pathway to produce SA from isochorismate in Arabidopsis, which is distinct from how SA is biosynthesized in bacteria. This study closes a major knowledge gap in plant SA metabolism and would help develop new strategies for engineering disease resistance in crop plants.

Independent evolution of rosmarinic acid biosynthesis in two sister families under the Lamiids clade of flowering plants.

As a means to maintain their sessile lifestyle amid challenging environments, plants produce an enormous diversity of compounds as chemical defenses against biotic and abiotic insults. The underpinning metabolic pathways that support the biosynthesis of these specialized chemicals in divergent plant species provide a rich arena for understanding the molecular evolution of complex metabolic traits. Rosmarinic acid (RA) is a phenolic natural product first discovered in plants of the mint family (Lamiaceae) and is recognized for its wide range of medicinal properties and potential applications in human dietary and medical interventions. Interestingly, the RA chemotype is present sporadically in multiple taxa of flowering plants as well as some hornworts and ferns, prompting the question whether its biosynthesis arose independently across different lineages. Here we report the elucidation of the RA biosynthetic pathway in Phacelia campanularia (desert bells). This species represents the borage family (Boraginaceae), an RA-producing family closely related to the Lamiaceae within the Lamiids clade. Using a multi-omics approach in combination with functional characterization of candidate genes both in vitro and in vivo, we found that RA biosynthesis in P. campanularia involves specific activities of a BAHD acyltransferase and two cytochrome P450 hydroxylases. Further phylogenetic and comparative structure-function analyses of the P. campanularia RA biosynthetic enzymes clearly indicate that RA biosynthesis has evolved independently at least twice in the Lamiids, an exemplary case of chemotypic convergence through disparate evolutionary trajectories.

Sequenceserver: A Modern Graphical User Interface for Custom BLAST Databases.

Comparing newly obtained and previously known nucleotide and amino-acid sequences underpins modern biological research. BLAST is a well-established tool for such comparisons but is challenging to use on new data sets. We combined a user-centric design philosophy with sustainable software development approaches to create Sequenceserver, a tool for running BLAST and visually inspecting BLAST results for biological interpretation. Sequenceserver uses simple algorithms to prevent potential analysis errors and provides flexible text-based and visual outputs to support researcher productivity. Our software can be rapidly installed for use by individuals or on shared servers.

The biosynthetic origin of psychoactive kavalactones in kava.

Kava (Piper methysticum) is an ethnomedicinal shrub native to the Polynesian islands with well-established anxiolytic and analgesic properties. Its main psychoactive principles, kavalactones, form a unique class of polyketides that interact with the human central nervous system through mechanisms distinct from those of conventional psychiatric drugs. However, an unknown biosynthetic machinery and difficulty in chemical synthesis hinder the therapeutic use of kavalactones. In addition, kava also produces flavokavains, which are chalconoids with anticancer properties structurally related to kavalactones. Here, we report de novo elucidation of the key enzymes of the kavalactone and flavokavain biosynthetic network. We present the structural basis for the evolutionary development of a pair of paralogous styrylpyrone synthases that establish the kavalactone scaffold and the catalytic mechanism of a regio- and stereo-specific kavalactone reductase that produces a subset of chiral kavalactones. We further demonstrate the feasibility of engineering styrylpyrone production in heterologous hosts, thus opening a way to develop kavalactone-based non-addictive psychiatric therapeutics through synthetic biology.

The Reference Genome Sequence of Scutellaria baicalensis Provides Insights into the Evolution of Wogonin Biosynthesis.

Scutellaria baicalensis Georgi is important in Chinese traditional medicine where preparations of dried roots, "Huang Qin," are used for liver and lung complaints and as complementary cancer treatments. We report a high-quality reference genome sequence for S. baicalensis where 93% of the 408.14-Mb genome has been assembled into nine pseudochromosomes with a super-N50 of 33.2 Mb. Comparison of this sequence with those of closely related species in the order Lamiales, Sesamum indicum and Salvia splendens, revealed that a specialized metabolic pathway for the synthesis of 4'-deoxyflavone bioactives evolved in the genus Scutellaria. We found that the gene encoding a specific cinnamate coenzyme A ligase likely obtained its new function following recent mutations, and that four genes encoding enzymes in the 4'-deoxyflavone pathway are present as tandem repeats in the genome of S. baicalensis. Further analyses revealed that gene duplications, segmental duplication, gene amplification, and point mutations coupled to gene neo- and subfunctionalizations were involved in the evolution of 4'-deoxyflavone synthesis in the genus Scutellaria. Our study not only provides significant insight into the evolution of specific flavone biosynthetic pathways in the mint family, Lamiaceae, but also will facilitate the development of tools for enhancing bioactive productivity by metabolic engineering in microbes or by molecular breeding in plants. The reference genome of S. baicalensis is also useful for improving the genome assemblies for other members of the mint family and offers an important foundation for decoding the synthetic pathways of bioactive compounds in medicinal plants.

Monoamine Biosynthesis via a Noncanonical Calcium-Activatable Aromatic Amino Acid Decarboxylase in Psilocybin Mushroom.

Aromatic l-amino acid decarboxylases (AAADs) are a phylogenetically diverse group of enzymes responsible for the decarboxylation of aromatic amino acid substrates into their corresponding aromatic arylalkylamines. AAADs have been extensively studied in mammals and plants as they catalyze the first step in the production of neurotransmitters and bioactive phytochemicals, respectively. Unlike mammals and plants, the hallucinogenic psilocybin mushroom Psilocybe cubensis reportedly employs an unrelated phosphatidylserine-decarboxylase-like enzyme to catalyze l-tryptophan decarboxylation, the first step in psilocybin biosynthesis. To explore the origin of this chemistry in psilocybin mushroom, we generated the first de novo transcriptomes of P. cubensis and investigated several putative l-tryptophan-decarboxylase-like enzymes. We report the biochemical characterization of a noncanonical AAAD from P. cubensis ( PcncAAAD) that exhibits substrate permissiveness toward l-phenylalanine, l-tyrosine, and l-tryptophan, as well as chloro-tryptophan derivatives. The crystal structure of PcncAAAD revealed the presence of a unique C-terminal appendage domain featuring a novel double-ß-barrel fold. This domain is required for PcncAAAD activity and regulates catalytic rate and thermal stability through calcium binding. PcncAAAD likely plays a role in psilocybin production in P. cubensis and offers a new tool for metabolic engineering of aromatic-amino-acid-derived natural products.

Contribution of Untargeted Metabolomics for Future Assessment of Biotech Crops.

The nutritional value and safety of food crops are ultimately determined by their chemical composition. Recent developments in the field of metabolomics have made it possible to characterize the metabolic profile of crops in a comprehensive and high-throughput manner. Here, we propose that state-of-the-art untargeted metabolomics technology should be leveraged for safety assessment of new crop products. We suggest generally applicable experimental design principles that facilitate the efficient and rigorous identification of both intended and unintended metabolic alterations associated with a newly engineered trait. Our proposition could contribute to increased transparency of the safety assessment process for new biotech crops.

Genetic regulation of mitotic competence in G0 quiescent cells.

Quiescent (G0 phase) cells must maintain mitotic competence (MC) to restart the cell cycle. This is essential for reproduction in unicellular organisms and also for development and cell replacement in higher organisms. Recently, suppression of MC has gained attention as a possible therapeutic strategy for cancer. Using a Schizosaccharomyces pombe deletion-mutant library, we identified 85 genes required to maintain MC during the G0 phase induced by nitrogen deprivation. G0 cells must recycle proteins and RNA, governed by anabolism, catabolism, transport, and availability of small molecules such as antioxidants. Protein phosphatases are also essential to maintain MC. In particular, Nem1-Spo7 protects the nucleus from autophagy by regulating Ned1, a lipin. These genes, designated GZE (G-Zero Essential) genes, reveal the landscape of genetic regulation of MC.